Pseudorapidity Distribution of Charged Particles in pp Collisions at y/s = 630 GeV 
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Using a silicon vertex detector, we measure the charged particle pseudorapidity distribution over 
the range 1.5 to 5.5 using data collected from pp collisions at ^/i = 630 GeV. With a data sample of 
' 3 x 10 6 events, we deduce a result with an overall normalization uncertainty of 5%, and typical bin to 

Q\ , bin errors of a few percent. We compare our result to the measurement of UA5, and the distribution 

generated by the Lund Monte Carlo with default settings. This is only the second measurement at 
this level of precision, and only the second measurement for pseudorapidity greater than 3. 
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The mean angular density of particles produced in high energy hadron-hadron collisions is a fundamental property 
of the strong interaction. Multi-particle production is inherently a non-perturbative process, it is currently best 
treated theoretically using Monte Carlo models implemented in computer programs such as Herwig [El, Isajet 
and PYTHIA Besides its fundamental importance, the angular particle density is a primary consideration in 
experiment design. Estimates of trigger rates and occupancies for future LHC detectors depend on extrapolations of 
our knowledge of particle production from energies of less than 2 TeV to 14 TeV. Future experiments to search for the 
quark-gluon plasma will require an accurate knowledge of the pseudorapidity distribution of single nucleon collisions 
CT) | to compare with heavy ion results. 

We report on the angular density of charged particles produced in 630 GeV pp collisions at the CERN SppS collider. 
A convenient way to characterize the angular density is through the pseudorapidity distribution, dN c h/dr], where N c h 
is the number of charged primary particles per collision, and rj = — lntan(#/2), 8 being the polar angle from the 
^ ' collision axis in the center of mass. Previous measurements of dN c h/drj in this energy range have been reported by 
1} , the UA5 [0, UAl ||, and CDF || experiments. Of these, only the UA5 result is of high precision. It is also the only 
O 1 measurement covering | rj |> 3. 

The data reported here are from the test of a forward geometry silicon micro- vertex detector proposed as part of a 
hadronic B-physics experiment (P238) 0. The geometry of the detector and the trigger play an important role in this 
analysis. As shown in Fig. Fj, the detector consists of 6 tracking planes. Each plane is composed of 4 quadrants, and 
each quadrant contains 2 silicon micro-strip detectors (MSD's) oriented orthogonally. One measures the x coordinate 
and the other measures the y coordinate. The coordinate system of the experiment is oriented with z along the 
direction of the proton beam, x horizontal, and y vertical. The MSD's have a diode pitch of 25 /im with every second 
diode readout through a capacitively coupled electrode. The MSD's have an active area of 4.48 cm x 4.48 cm. 

The overall length of the detector is less than 20 cm. Due to its small size relative to the longitudinal extent of the 
luminous region (er ps 12 cm), we reconstruct particles with 1 < \r)\ < 6, depending on the longitudinal position of the 
primary interaction (see Fig. ||). However, the detector covers about 3 units of r] for any given position of the primary 
interaction. The positions of the MSD's are determined from analysis of the data to a precision of a few microns in 
x and y and about 100 //m in z. Q First, the MSD's within a quadrant are aligned by optimizing track residuals. 
Second, the quadrants are aligned to each other by optimizing vertex residuals. In the second step, the z-axis of the 
coordinate system is forced to coincide with the axis of the luminous region. The procedure is described in detail in 
Ref. [||. After alignment, we are able to determine the spatial distribution of primary interactions to high accuaracy. 

About 5 million events were recorded to tape with a minimum bias trigger. The trigger consists of 2 scintillation 
counters placed roughly ±3 m from the detector. A coincidence of the 2 counters in time with the beam crossing 
forms the trigger. 

Event reconstruction begins by forming clusters from adjacent fired strips. Tracks are formed by grouping clusters 
on different planes that lie within a road of 50 /xm half width fjj. The x-z and y-z track projections (referred to as 
x view and y view tracks) are reconstructed independently. The two views are coupled for the purposes of finding a 
vertex, since the z coordinate of the vertex is common to all the track projections in an event. Vertex finding begins 
by determining the intercepts of the track projections with the beamline. A cluster of intercepts from at least 2 x and 
2 y view tracks within 1 cm forms the seed for a vertex. We assign all tracks passing within 600 /im of the position of 
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the seed to the vertex. The process is repeated to search for other vertices using the tracks not previously assigned 
to a vertex. 

We derive dN c h/dr) from the data in two ways: from the x and y view tracks, and from a subset of track projections 
that are matched in x and y to yield a 3-dimensional track ||. Only events with at least one reconstructed vertex 
are selected for inclusion in the dN c h/ ' dr\ analysis. Only tracks assigned to a vertex are used in the analysis, removing 
most spurious tracks as well as many tracks from secondary decays (primarily K®'s and A's). 

Letting u stand for x or y, we define the r\ projections r/ u = — lntan(# u /2), where 9 U is the projected angle of 
the track with respect to the beam axis. Using u view tracks, we form a binned dN c h/drj u distribution whose i 
bin covers the interval r\ % u to rf u + Aij u . The raw dN c h/dn u distribution from the data is corrected, using the Monte 
Carlo simulation, for any remaining spurious or secondary decay tracks and tracks from single diffractive interactions. 
The corrections for spurious tracks and tracks from secondary decays are each about 2%. The single diffractive 
interaction correction is about 0.5%. These corrections have roughly the same shape over the measurement interval 
as the corresponding dN c hjdr\ u distribution. 

The simulation is used to unfold dN c h/dr) from dN c h/dr) u by minimizing 

X 2 = (A u e - gf V" 1 (A u e- g) , (1) 

where e represents the unnormalizcd result vector for dN c h/drj, and g represents the unnormalized vector of measure- 
ments. The matrix A„ is determined from the simulation and accounts for the effects of acceptance and smearing as 
well as projecting dN c h/drj. The (ij) th element of A„ is the fraction of charged primary particles generated with ij 
in bin j that are reconstructed and have r\ u lying in bin i. Simulated events are reconstructed with the same code as 
actual data. The covariance matrix, V ff , accounts for the finite statistics of the data and simulation. 

The second method determines dN c h/dr) directly (without unfolding) using tracks reconstructed in 3-dimensions. 
We match x and y view tracks whenever they traverse a unique combination of MSD's. The probability that a match 
can be made depends on the overall distribution of tracks in the event and is determined by simulation. An additional 
source of background arises from incorrect matchings. We simulate and correct for this background as well as spurious 
and secondary tracks. Eq. [I] is used to determine dN c h/dr) from the matched tracks, where A contains the effects of 
acceptance, smearing, and the probability of matching. 

The quantity e which minimizes Eq. |l]is proportional to dN c h/dr], 

dN ch _ N MC g ^ 
drj N gen N data Ari 

where N^ata is the number of vertices reconstructed in the data, Nmc is the number reconstructed in the simulation, 
N gen is the number of generated interactions, and An — 0.25 is the bin width. 

A good simulation is needed to make the corrections for acceptance and reconstruction efficiency. We use a detailed 
detector simulation based on GEANT 3.15, coupled to PYTHIA 5.7 and JETSET 7.4 for event generation. All the 
standard features of GEANT are used: delta-ray production, Moliere model of Coulomb scattering, decays of long 
lived particles, photon conversions, hadronic interactions, etc. The simulation of the MSD's is done at the individual 
diode level, and includes charge sharing between diodes and detector noise. The positions of the MSD's used in the 
simulation are determined from the data. During analysis, signals on non-functioning channels are suppressed, and 
random detector inefficiencies are introduced, according to the detector performance for the data set being simulated. 

Events are generated with a spatial distribution of primary interactions and pile-up rate determined from the 
data p0| . The charged particle multiplicity distribution is generated according to the measurement of UA5 [Q. 
To simulate the UA5 distribution, we first enable multilple interactions in PYTHIA, then weight the events. The 
scintillation counters comprising the trigger and the section of beampipe between the vertex detector and the counters 
are included in the simulation . The trigger is simulated by requiring a minimum energy deposition in each counter 
for the event to be triggered. Due to the trigger, the probability for finding a track has a maximum around \n\ of 3. 

The results for approximately 3 x 10 6 triggers are shown in Fig. [}| The results are shown for x view, y view, and 
matched tracks. The agreement is good between the x and y view points up to \rj\ — 4. Furthermore, the shapes of all 
3 measurements are quite similar. The differences in the shapes are used to determine the errors, as described below. 
The results of Fig. § remain consistent when subdivided by direction of track (±z), detector quadrant, or number of 
detector planes each track traverses . 

We calculate the final result as the average of the unmatched (x-view and y-view), and matched track distributions. 
The error has 2 contributions: an overall normalization error of 5% (half of the 10% difference between matched and 
unmatched), and individual bin errors relating to the shape of the distribution. While the normalization error is 
entirely systematic, the shape error has both systematic and statistical components. The shape error is determined 
from dN c h/drj results obtained from 6 independent subsets of the data. The subsets consist of data collected during 
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particular running periods, and contain comparable statistics. The average of the matched and unmatched results is 
found for each subset. The deviations of these averages from the final result are displayed in Fig. [|. Treating each 
as statistically equivalent, the shape error is estimated from the RMS spread of the points. Note that the excellent 
agreement of these results is obtained with no adjustments to the normalization, attesting to the consistency of the 
result. Half the difference between the x and y results is added in quadrature to obtain the total shape error. 

We find a number of other possible systematic errors to be negligible in comparison to the systematic errors cited 
above jfl 12|. These include uncertainties due to the detector efficiencies, the locations of the MSD's, the positions of 
the trigger scintillators, the charged multiplicity distribution, the dN c h/drj distribution generated in the simulation, 
the single diffraction cross section |u|, and the requirements for matched tracks. We also calculated dN c h/dr] using 
subsets of the data which might be sensitive to systematic errors. The subsets consist of: tracks intersecting a 
specific number of planes, sensitive to errors in detector efficiencies and misalignments; tracks in a particular detector 
quadrant, sensitive to the relative alignment of MSD's within a quadrant; and tracks with positive (negative) values of 
77, sensitive to the longitudinal distribution of interactions. From these data subsets we find no indication of additional 
systematic effects. 

The dN c h/drj values and the shape errors are listed in Table Q. Fig. || displays our result along with the UA5 
result Q and the distribution generated by PYTHIA, unweighted for the purposes of comparison. The prediction 
of PYTHIA agrees quite well with our result in shape, but is low in normalization. Our agreement with the UA5 
measurement is reasonable, however there is an indication that our measurement falls off less rapidly at high \rj\. 

In summary, we have made the second high precision measurement of dN c h/drj at a high energy hadron collider. 
The most significant systematic error is the 5% normalization error arising from the difference between the matched 
and unmatched track results. Errors related to the shape of the distribution are at the level of a few percent. We find 
good agreement with the shape prediction of the Lund Monte Carlo. We find reasonable agreement with previous 
measurements, though there is an indication that our measurement falls off less rapidly at high \rj\ compared with 
UA5. 
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TABLE I. The values of dN ch /dri and their shape errors. The 5% normalization error is not included. 



r\ range 


dN ch 1 dr\ 


error 


1.50-1.75 


3.30 


0.16 


1.75-2.00 


3.23 


0.11 


2.00-2.25 


3.30 


0.07 


2.25-2.50 


3.25 


0.04 


2.50-2.75 


3.16 


0.04 


2.75-3.00 


3.02 


0.03 


3.00-3.25 


2.88 


0.02 


3.25-3.50 


2.71 


0.03 


3.50-3.75 


2.57 


0.04 


3.75-4.00 


2.40 


0.05 


4.00-4.25 


2.24 


0.09 


4.25-4.50 


2.07 


0.12 


4.50-4.75 


1.84 


0.16 


4.75-5.00 


1.68 


0.20 


5.00-5.25 


1.45 


0.32 


5.25-5.50 


1.19 


0.47 
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FIG. 1. Isometric view of the layout of the MSD's of the P238 test. The p and p beams pass through the gap between the 
upper and lower sets of MSD's. The horizontal and vertical lines represent the directions of the strips on the MSD's. Each 
MSD is 4.5 x 4.5 cm and planes are separated by 3.8 cm. 
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FIG. 2. Side view of the upper half of the detector showing the acceptance limits for tracks from vertices at z — —45 cm and 
— 10 cm. The luminous region is centered around z — —8 cm with an rms of about 13 cm in z. For this analysis, we use events 
produced within —45 cm < z < 30 cm. Although the acceptance limits change by about 1.5 units of t] between the 2 positions, 
the total coverage is about 3 units of r\ for each. 
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FIG. 3. dN c h/dri derived from dN c h/drj x , dN c h/dr]y, and dN c h/ 1 dr] m . The bin width is 0.25 units of r\. Only statistical errors 
are shown. 
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FIG. 4. The differences between dN c h/dri calculated separately for 6 data subsets and the final result. The differences arise 
from a combination of statistical and systematic errors and represent part of the shape error quoted in the final result. 
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FIG. 5. Final dN c h/dri result compared to the result from UA5 and the distribution generated by PYTHIA-JETSET. For 
the P238 points, only the shape errors are shown, the 5% normalization error is not plotted. The UA5 points are plotted 
with their statistical errors, a normalization uncertainty of order 2.5% is not shown. The Monte Carlo curve is unweighted for 
purposes of comparison. 
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